m ■ 

o : 

A new scheme for the Klein-Gordon and Dirac 
o : fields on the lattice with axial anomaly 

<D ■ 

^ ■ Miguel Lorente 

. Departamento de Fisica, Universidad de Oviedo, 33007 Oviedo, Spain 

Arnold Sommerfeld Institut fur Theoretische Physik, Technische Universitat, D-3392 Clausthal, 
l ' Germany 

. (Submitted 30 October 1992) 

o : 

Abstract 

Using the method of finite differences a scheme is proposed to solve exactly the Klein-Gordon 
and Dirac free field equations, in a (1 + l)-dimensional lattice. The hamiltonian of the Dirac 
field is translational invariant, hermitian, avoids fermion doubling, and, for the massless case, 
»- H , preserves global chiral symmetry. Coupling the fermion field to the electromagnetic vector 

I potential we construct a gauge invariant vector current leading to the correct axial anomaly. 

D ! PACS number: 02.20.+b, ll.30.-j 
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;rj ■ 1 Introduction 

rS ' 

The method of finite elements has become very powerful to solve evolution equations in quantum 
field theories. In particular Bender and his collaborators (l| have applied the method to the 
Klein-Gordon and Dirac quantum field equations, consistent with the equal time commutation 
relations, which preserves global chiral symmetries in the massless case and avoids fermion doubling. 
Matsuyama [Q] has constructed an abelian gauge-coupled Dirac equation without species doubling. 
Milton H extended this method to the non-abelian gauge theories. Using the method of finite 
differences, Vazquez j|] also introduced an explicit scheme for the Dirac fields, such that the 
associated action is hermitian, the chiral symmetry is preserved and the fermion doubling is absent, 
but the equal time commutation relations are not satisfied. Nevertheless in all these models the 
solution is given for only one time step in the transfer matrix. 

Lattice field theories have become very popular to analyze perturbation theories in the standard 
model H . With the lattice regularization comes the problem of fermion doubling if one imposes to 
the hamiltonian translational invariance, hermiticity and locality. Under these assumptions fermion 
doubling is unavoidable, as Nielsen and Ninomiya have proved ||. 

We present a simple model for the Klein-Gordon and Dirac free field using the method of finite 
differences. The novelty of our approach lies, first, on the exact solutions for any time step of the 
difference equations, making a good contact with the covariant continuous theory (we have already 
worked out exact solutions for the Heisenberg difference equations of motions in a time lattice J?]). 



Secondly, we try to escape the no-go theorem of Nielsen and Ninomiya because our Hamiltonian is 
non-local, although avoids fermion doubling. 

To check our model we apply it to current algebra, namely, we implement the vector and 
axial current of fermion massless fields with abelian gauge fields. The axial current is then gauge 
invariant, but exhibit non-local point separation structure, which in the continuous limit leads to 
the correct axial anomaly. 



2 A set of orthonormal functions on the lattice 

Let us define the function 

(l + UkeV 

where k is an arbitrary constant and e some small quantity. 
This function satisfies the difference equations 

Sj+i - Si = ike \ Ui+l + Sj) , fj ~ fj-i = ike \ iSi + fj-i) ( 2 ) 

Using the notation for the difference operators 

A fj = Si+l ~ fj , V/j = fj - fj-i , Sfj = f j+ i - f._i 
and the average operators 

~ A Si = \{fi+i + Si) , ^Si = \ (Sj + Si-i) , nS3 = \(s j+ % + Si-\) 

the above equations become 

jAfj = ikAfj^Vf^ifiVfj (4) 
-Sfj = ikfifj (5) 
These difference equations and their solutions are approximation of the differencial equation 

(6) 

with truncation error of second order in e (the truncation error is the difference between the Taylor 
expansion of the finite difference equation and the differential equation). 

On the other hand the function fj converges to the solution of the differential equation (|6|) : 

when j — » oo, je — > x, e — > with order of convergence O (e 2 ) . 



If we impose the boundary conditions /o = /at = 1 for some fixed number N, we find N different 
values of k 

2 irm , . 

/c m = -tan— , m = 0,l,...,7V-l (8) 



with properties k m +N = k m and k. 

fj [km) = I 



m = —k m - Therefore 

\ j 

m = 0,1,..., N- 1 



1 -\~ <2^Skyyi\ 



1 ^tskffi 



(9) 



Notice that if m = N/2 for TV even, then fj [k m ) = (—1) J ■ Obviously, 



fj (k m ) = exp i—mj 



2vr 



N 



These functions satisfy the orthogonality relations 



N-l 



~pr fj (km) fj (k m i) — 5 n 



(10) 



(11) 



3=0 



with respect to the scalar product in the Hilbert space <jL ([0, TV]). With the help of this orthog- 
onality relation we can express any function Fj of discrete variable in terms of the TV" orthogonal 
functions fj (k m ) , namely, 



N-l 



^ ' Q"mfj (kr, 



m=0 



with 



N-l 



"Trf E # M F i> m = 0,l,...,N-l 



The orthonormal set 



fj (km, 



N-l 



m=0 



(12) 



(13) 



(14) 



can be chosen as a basis for the Hilbert space (jL ([0,7V]) , hence, a completness relation can be 
constructed 



1 ^ 1 



TVe 



y \ n m J 



S; 



33 



(15) 



m=0 



the right hand side being the discrete analog of the Dirac Delta function. 
Consider now the following function on the lattice 



«* = 3^Er {/i ^) - 1} 

m=l " l 



(16) 



s=N 



Using the difference equations @ we can calculate Uj, namely, 

2 E 7^{/i-/i-l + / J -l-/ J -2 + ... + /2-/l + /l-/o} 



TVei ^ fc 

m=l 



2 E^ £ U, + /^l + / J -l + /,-2 + -.. + / 2 + /l + /l + /0} 



TVei ^ 2fc 

m=l 



1, Vj > 



1 


j > 




1/2 


3 = 


(17) 





j<0 





Similarly, Uj = — 1, Vj < 

Hence we can have the discrete analog of the Heaviside function, 

#3 = 2 («i + X ) 
It is easy to prove the difference equations 

\ Mj = Ye^ ]+lfi + 5jfi)l 7 W i = 2^ + ^-i,o) (18) 
\ M * : 2 £ 



ifeo + ^-io) (19) 



conecting the discrete versions of the Heaviside and Dirac functions. 
We have also the useful relations: 



A(/ i5i ) = (Afj) (A 9j ) + (A/,-) (A 9j ) (20) 
A(f j9j ) = (A/,) (A 9j ) + 1 {Afj) (A 9j ) (21) 

3 Hamiltonian formalism of the Klein-Gordon field on a (1 + 1)- 
dimensional Minkowski lattice 

We introduce the method of finite differences for the Klein-Gordon scalar field. An explicit scheme 
for the wave equation consistent with the continuous case (the truncation error is of second orden 
with respect to space and time variables) can be constructed as follows: 

^VnAnVjAj - ^VjAjVnAn + M 2 V n A n V j Aj^j <P] = (22) 

where the field is defined in the grid points of the (1 + l)-dimensional lattice 0™ = </> (je, nr) , e, r 
being the space and time fundamental intervals, j,n integer numbers and Aj (Vj) are the forward 
(backward) differences with respect to the space index, Aj (VjJ the forward (backward) averages, 
and similarly for the time index. 

Using the method of separation of variables it can easily be proved that the following functions 
of discrete variables are solutions of the wave equation (22): 

provided the "dispersion relation" is satisfied: 

to 2 — k 2 = M 2 (24) 

M, being the mass of the particle. 
We have also the solutions of (p2]) 

</>] = ("l) J+n (25) 



In the limit, j — * do, n — > do, je — > x, nr — > t the functions (|23| ) become plane wave 
solutions 

ff (k, uj) -» exp i (kx - cot) (26) 



Imposing boundary conditions on the space indices 
we get 



therefore 



2 Trm 

fc m = -tan— , m = 0, 1, . . . ,N - 1 



1/2 



For the positive energy solutions we define 



(27) 

(28) 
(29) 

(30) 



The solutions fl25l) £11*6 ct particular case of /j 71 (k m ,iv m ) with m = iV/2. 

Starting from the wave equation (^2[) we can construct a current vector. Multiplying ([22] ) 
by V n AnV jAj(J)j U from the left, and multiplying the complex conjugate of the wave equation by 
V n A n VjAj</>" from the right, substracting both results and using (pp|) we obtain the "conservation 
law" 

-Vjii - i-Vnk = (31) 



where 



Ji 

34 



- £ A 3 - (v n A n <f>* n ) Aj (v„A n <#) - Aj (v n A n $ n ) ~Aj (v n A n <^ 
-A n (VjA^f) A n (VjArff) - A n (V.Ajuf) -A n (V.A.o'j 



(32) 
(33) 



are the spatial and time component, respectively, of the charge vector current on the lattice. 

The charge density p suggest that we can substitute the scalar field 4>(x,t) by the promediated 
quantities VjAjft- and (jf (x,t) by VjAj(f)* n . 

Now we address ourselves to the real Klein-Gordon field . A suitable Hamiltonian for the field 
(fr™ and its conjugate momentum 7r™ can be defined as follows 



N-l , , 
3=0 * k 



rt) +A v ^^l ] +M 



(34) 



which is obviously hermitian. Using the periodicity condition of the fields and the identity (^) we 
can integrate by parts the second term on the right with the result: 



N-l 



j=0 



N-l 



Hn = eY, \ {iyA^lf ~ 4 ( V / A /"/) (?A^) + M 2 (V 3 A^) 2 \ (35) 



3=0 



As in the continuous case, we can derived the Hamilton equations of motions, varying the 
Hamiltonian density TCj first with respecto to the promediate momentum and secondly with respect 
to scalar field: 



1 

r 



-A n IVjArf] 



A n (vfArf 



dH n f 1 



n 
3^j 



(36) 

(37) 



(The same result was obtained by Bender and Milton (1986) if one eliminates the auxiliary field 
r in their formule (2.5) of Reference 8). 

Applying the difference operator ^V n on both sides of ( |3tf ) and substituying ([37]) in the result 
we recover the wave equation (|22]). 

Using (|3^) and (|37| ) it can easily be proved that the Hamiltonian ( |35|) is independent of the 
time index re, namely: 

V n H n = A n H n = (38) 

hence 

V n i? n = A n H n = H n (39) 

Since the plane wave solutions /j 1 (k m ,uj m ) (m = 0, 1 . . . N — 1) form a complete set of orthog- 
onal functions, we can expand the wave field and its conjugate momentum as 



m=—N/2 v 



71"; 



We 



E ( X + ^^m) ( a m/™ {km, U m ) - a* m f* n (k m ,U m ) 



N/2-1 



m=-N/2 

Applying the average operator on both sides we get 

N/2-1 



(40) 
(41) 



Ne m=-N/2 ^™ V 

_. N/2-1 r— 

\l ~2~ (° m /j (^mjk-Vn) ~~ a mfj {krni^rn) 



'Ne 



(42) 



(43) 



m=-N/2 

Using the completness relation (FL5|) we can calculate the coefficients in the Fourier expansion 



N-i 

E 



iV-1 



^f* n (k m ,u; m ) VjArf? + ^=/; n (k m ,u; m ) VjArf 

v 



we u 



„_n I V ^ V ZuJ m 



(44) 
(45) 



Using (|20D, the periodicity condition and Hamilton equations of motions (36, 37), one can easily 
prove that this coefficients are independent of re, namely, 

V n (o m ) = A n (am) = V n («4) = A n (a* m ) = (46) 

Turning to the complex Klein-Gordon field we can define a Hamiltonian as follows: 

N-l , 1 v 

Hn = e^2\ (ViA i7 r*) - (v.A^f ) V,A^ + M 2 (v, A^f ) V,A,^ (47) 

3=0 <- b ' 

from which the Hamilton equations of motion can be derived: 

-A„\ A/') = AnfijAjTT*/ 



-A n [VjAjTv] 



An [ pVjA^f 



M 2 (V, A 



(48) 
(49) 



leading to the wave equation (p2[). The solution of the wave equation can be expanded in Fourier 
series as before: 



! n 
3 



1 N ^ (l + \e 2 kl 



7T; 



'Ne 



HVe 



^2 



Q"mfj fomi^ra) ^mfj former. 



E 

m=-N/2 
N/2-1 

E 1 + 4 A ™ (b m f?(k m ,LO m )- a * m f™(k m ,u m )) 

m=-N/2 V V ' 



and similar expressions for the complex fields. 

In order to make connection of our scheme with the Einstein-de Broglie relations E = Huj, p 



hk we take for the period T and wave length A of the discrete plane waves functions ( 23 ) and (|27| ) 

T = Nt, A = Ne 

and for the phase velocity 



A e 
v p = T = t 

We have define the wave number and the angular frequency of the wave functions as: 

2 7rm 2 -Km 

k m = — tan — — , uj m = — tan — — , m = 0, 1, . . . , N — 1 



2 7rm 
-tan—-, o;„ 

e iV 



2 7rm 
-tan—, 
r W 



substituting the Einstein-de Broglie relations in the relativistic expresion E 2 — p 2 = M 2 (we use 
natural units h = c = 1), we obtain 



2 — h 2 — 2 



1 



M 



Since the phase velocity and group velocity satisfy v p v g = 1, we have finally 

M 2 



4 Quantitation of the Klein-Gordon field 



Let us start form the Hamilton equation of motion (36, 37) where the real fields and 7r™ now 
are, in general, non commuting operators: 



n 
3^3 



A n (VjArf 
A n (pV^ - M 2 [vfA^f] 



(50) 
(51) 



Eliminating the field 7r™ we obtain again the wave equation on the lattice (22) the solutions of 
which can be expanded in terms of the "plane wave functions" 



N/2-1 



j-r- E/ /o (°m/j {kmi^m) + a mfj for. 



'Ne \/2u. 

m=—N/2 
N/2-1 

—I 



'Ne 



m=-N/2 



, Ul r , 



(52) 
(53) 



where the coefficients a m , a m , are, in general, non-commuting operators. 

We introduce the quantization procedure imposing the commutation relations of the fields for 
n = 



(54) 
(55) 

(56) 



In order to prove unitarity we must prove that these relations hold also for any n. Since (|5 
and (1537) are valid for n = 0, we solve them for a m and 



1 \0 
\CL m ) 

(at) 



N-l , 1 



'Ne 
e 

We 



J 

N-l 



j=0 



V2 



Using (p4]-56) it is straight forward to prove that 



( a l)°>( a L 



(am)°, 

From ([5(], 51) and (57, 58) one can prove (remember [!(]) 
hence (59, 60) is also true for n = 1, and from ([521 |53|) we derive 



•a 



and the proof is completed by induction. 

The Hamiltonian for the real scalar field reads: 

N-l 



H n = * E \ {(VA-^) 2 - 1 (v.A^) (V,A,o'/) + ^(VA^ 



j=0 

A convenient scheme for the Heisenberg equation is: 



-A n (Vj-A^i 



1 

T 



-A n (VjAjt" 



A„ V,-A 



A„ |V,-A,-7r? 



Using (^l|) and ( |39| ) and ( |50| , 51) for any n, we get for the first equation 



-A n (V ; A ; o;') = -A n [VjAjfyHn] = A n (VjAj-irJ 



(57) 
(58) 



(59) 
(60) 

(61) 

(62) 
(63) 



(64) 



(65) 



and similarly for the second equation 



1 



-A n IVjArf) = A n -^VjA^n - M'VjAjfl 



(66) 



consistent with the Hamilton equation of motions (50, 51). 



Substituting (]57], |5|) in (||) we find 



N/2-1 



CL /m Q>r 



m=-N/2 

If we define the linear momentum on the lattice as 

N-l ^ 



(67) 



p = c E - PM) 



j=0 



we find 



N/2-1 

P — ^ ' k m {^a, m (i^ ri -\- ci^cLf 

m=-N/2 



N/2-1 



5 ' km^m^r 



m=-N/2 

because the zero point momentum vanish by cancellation of k m with. k— m 
For the complex Klein-Gordon field we have the hermitian Hamiltonian 

H n = e ^{(ViA^t) (v^) - p (ViA,-^) +M 2 (v i A i 0t n ) (v^)} 

(68) 

(69) 



j=0 



The equal time commutation relations now read: 



-din' 



33 ' 



with all other equal time commutators vanishing. 

The Heisenberg equations of motion lead, as in the real case, to the discrete wave equation for 
the fields 7r" , and their adjoints. The Fourier expansion in terms of the plane wave solutions 
are 



1 N ^ (l+i---A-,„ 



m=—N/2 
N/2-1 



(a m fj (k m ,uj m ) + b\j* n (k m , u m f) 



7T; 



'Ne 



E (l + J^m) (k m , u m ) - atf* n (km, L0 r , 



m=-N/2 

Inverting this equations we can calculate with the help of (|69|). 



,t 



■"mm' i 



(70) 
(71) 

(72) 



with all other commutators vanishing. 

The Hamiltonian can be expressed in terms of these creation and annihilation operators: 



N/2-1 

H = 2 E W ™ ( a t«»n + 4 



(73) 



m=-N/2 



The charge of the complex field is given by 
N-l 



Q = -ieJ2 (VA'O (VA^) - ( VA-ttJ" ) ( VjAj< 

3=0 



N/2-1 

(a m a m + btb m ) (74) 

m=-N/2 



Finally, two very important quantities can be evaluated on the lattice: the general commutations 
rules and the Feynman propagator. 

First of all, using ( |70|) and ( |72|) in order to calculate the commutator of the fields for two 
different space-time points, we find: 



N/2-1 



(75) 



which becomes zero for equal time (n = n') 

Secondly, the Feynman propagator function is defined to be the vacuum expectation value 



A, 







(76) 



of the time ordered product 



n J n ' 



for n > n' 



for n' > n 



(77) 



Notice that there is no ambiguity of the time ordered product for equal times, sine and </>j 
commute for n = n' because of (f75|). 

The discrete analog of the Feynman propagator is from (|70|) and (|77|) 



fn' 



1 N/2-1 fl + i £ 2 k 



m=-N/2 



2iO r , 



{finfj 1 (k m ,LO m ) + &- n fj (k m , U m )\ 



(78) 



•d n being the discrete realization of the Heaviside functions jl7|). 

In order to see the analogy with the continuous case we could take the discrete Fourier transform 



of (|7q). Instead, we take the Fourier transform of the wave equation (22), namely, 

1 1 



(-<4 + k 2 m + M 2 



i + k 2 e i + \t 



1-2, ,2 Ym 







■a;, 



where 



N-l N-l 
j=0 n=0 

is the Fourier transform of the Klein-Gordon field. 

Therefore the Feynman propagator is defined in the momentum space 



and in the configuration space 



N/2-l N/2-l 



s-m—n 



(Ne) (Nt) 



E E & m fE?' 



m=-N/2r=-N/2 



(Notice that now u r = ^ tan ^£ , can take negative values) 

It can be proved that this propagator is the Green fuction for the discrete wave equation, 
namely, 

-j V n A„V j A i - -jVjAj-VnAn + M 2 V n A n VjAj\ G]Zj> = ~5 n n>-$jj> 
With the aid of the identities: 



V„A„ 



V n A n 



tfn/j 1 ] = tfnV n A n /; + - [(A n tf n ) (A n /;) + (A n tf n _x) (A^j 1 - 1 )] + - (V n A n tf n ) /" 

and the equations fli~8|) and (pjj) one can prove that Ap defined by (|78|) satisfies also the same 
condition, therefore Ap and are equal up to a constant. 



5 Quantization of the Dirac field 



The Hamiltonian for the Dirac field ipaj an d its hermitian adjoint ij)^ can be defined as: 

N-l 



H n = s E A^j" | 7 47i -A,V" + M 7 4A^™) 
3=0 1 £ J 



with 



71 



-i 

1 



74 



1 



, ^7l74 = 75 



-1 
-1 



(79) 



(80) 



In the next chapter it will be shown that the substitution tp (x) — > Ajip™ is consistent with the 
density function 

p = {h^Y) ( A ^" ) ^ 8i ) 

If we defined the conjugate momentum 7r™ = iifij the Hamilton equation of motion are obtained 

tn 

variying the Hamiltonian density H.™ with respect to iipj , that is to say: 

X -A n A^ n 3 = + M 74 A i ^| (82) 

from which it can be derived the Dirac equation on the lattice: 

Tl-AjAn - i7 4 -A n A j + MAjA n ) V" = (83) 



The same result was obtained by Bender, Milton and Sharp (Reference 1, formula 18) applying 
the method of finite elements to the action. 

From (EH) and (|83[) one proves for the Dirac Hamiltonian (|79|) 



A n H n = 0, hence A n H n = H n 



(84) 



If we impose periodicity conditions in the fields we can prove by "integration by parts" of 
(|79|)that H n is an hermitian function. 

Applying the operator 7i^VjV ri — ij^VnAj — MVjV n on both sides of ( |82] ) we recover the 
wave equation (]2^) for ip^ 

Let us construct solutions to ( |83| ) of the form 

W = w(k,E)f?(k,E) (85) 



where /j 1 (k, E) is given by fl23[). 

The four component spinors w(k,E) must satisfy 

(i^k -j 4 E + M)w(k,E) = (86) 

Multiplying this equation from the left by (i'fik — 74-E — M) we obtain the "dispersion relation" 
for the Dirac equation 

E 2 — k 2 = M 2 (87) 

Following standard procedure we can construct a complete set of Dirac wave functions of the 
form 



1 M 



— — 11 ,-, u mfj (km,E v 

'Ne V E m 



7=\ -^—Vmfj™ 1 (k m , E m ) 
VNe V E m 

where 

2 vr 



E m = +y k m + M 2 , k m = -tan—m, m = 0,±l,... (90) 



and the spinors u m and v m are defined as usual: 



1 



Um = w(km,E m ) = ( Em ^ M S ) [ km J ( 91 ) 



2M 



-Em + M \ 2 



W (-km, -E m ) = {^^) ' { Em l M j (92) 

With the aid of this orthonormal set we can expand the fields in the usual way (remenber f42| ) 

N/2-1 



A^" = -7= E \ -E^(umCmf?(km,E m )+V m dif* n (km,E m )) (93) 



m=-N/2 
N/2-1 



K 3 iT = -J= J^(umcif* n (k m ,Em)+Vmd m f?(km,E m )) (94) 

m=—N/2 ' 

where ^ = V , j™74- Inverting this expressions, we find 



7^\h^ £ £ /; n (&m,£m)«m74A^ (95) 
'JVe V , =0 



1 



M 



N-l 



fNe V K 



1 



M 



fNe V ^ 
1 



M 



/ iVe V 



£ kjlpjliUmf 1 - (k m , E m ) 

3=0 
' N-i 

£ Y &ji>jl4V m fj n (k m ,E m ) 

3=0 
' N-l 

£ E /™(^m,^ m )Wm74A j V" 
3=0 



(96) 

(97) 
(98) 



For the quantization of the Dirac field we start from the Hamilton equation of motion (|82| ) 

where now ij)^ and ipj™ are in general non commuting operators. Using the same technique as in 
the Klein-Gordon field we can prove that the equal time anticommutation relations hold for any n 

1 



(99) 
(100) 



provided it holds for n = 0. 

Introducing the Hamiltonian (^) in the Heisenberg equation of motion and using the anticom- 
mutation relations (99, |100| ) we recover the Dirac equation, namely, 



1 7 



A, 



Aj^,H n 



A n \ 7471- Aji/>2 + M lA A 3 ^ 



Using (99, 100) we can derived the anticommutation relations for the creation and annihilation 
operators 



t 



d (V 



(101) 
lation 

(102) 



with other anticommutation relations vanishing. 

Finally, the Hamiltonian (|79|) can be written in terms of these operators 



JV/2-1 

H — ^ ] E m (c^ n c m + d^dm 1 



(103) 



m=-N/2 

Similarly the momentum and charge operators are written: 

N-l , , x N/2-1 



j=0 
N-l 



£ (A^J n J (Arff) = Y k m (ctc m + dtd ri 



q = E(4-#)(4-^)= E (ctc m -did m ) + 



m=-N/2 
N/2-1 



2N 



The discrete analog of the Feynman propagator for the Dirac field is defined as before 



cm— n 

where the time order operator acts as follows 



i/)j"ipj> for n > n' 
—Wj' f° r n' > n 



(104) 
(105) 

(106) 

(107) 



Following standard methods we find: 

N/2-1 f± _ i_ £ fc \ 2 

S^'-f = ^We £ 2lE ^ k m + iE m + iM) af3 f^p'(k m ,E m ) 

m=—N/2 

N/2-1 h- lisJe Y 

- $-nJf £ E 2iE (-yi k m + iE m + iM) ap f^J, n '(k m ,E m ) (108) 

m=—N/2 m 

We can also derive the Feynman propagator from the Dirac equation. Taking the Fourier 
transform of (^) we get 

_ (i-jfefc.) (i + l„-E r ) _ , _ i , , i , (109) 

i7ifcm-74-Br + M V 2 "j \ 2 'J Ef - kl - M 2 

hence 

JV/2-1 JV/2-1 

^-7' = ^ E E frraftfi^E,.) (110) 
m=-N/2r=-N/2 

Both function ( |108| ) and ( |110| ) are the Green function for the Dirac equations, therefore they are 
equal up to a constant. 

Our model for the fermion field satisfies the following conditions: 

i) the hamiltonian ( f79[ ) is traslational invariant with respect to the space indices. 

ii) the hamiltonian is hermitian. 

iii) for M = 0, the wave equation (|83|) is invariant under global chiral transformations. 

iv) there is no "fermion doubling" as it can be seen in ([X]). In fact, E m takes the value M at 
m = and nowhere else. 



v) the hamiltonian is non-local. Using the finite Fourier transform ( |13| ) the hamiltonian H {krn) — 
i747ifc m +74M, given by (|86|), and consequently the dispersion relations (|9C|) are smooth func- 
tions of k m , except for m = N/2. Therefore our model escapes the no-go theorem of Nielsen 
and Ninomiya fj. 



6 Anomalies in axial vector current 



In order to construct vector and axial currents in terms of the fermions fields in the presence of 
external electromagnetic vector potential, we multiply ( p3| ) from the left by AjA n ipj and then we 
multiply the adjoint equation of (83) from the right by AjA n ^>™. Adding together both results and 
using (f20|) we find: 

±Aj (A n ^; 7l A n V™) - \iA n (A^^A^) = (111) 

This equation can be considered the discrete version of the "conservation law" for the vector 
current 

h = i (A n V^) 71 (A n ^ n ) , h = i (A^) 74 (A 3 ^) (112) 



The same equation (111) can be applied to the axial current 



,•5 



i (A n ^) 7175 (A n ^ n ) , jt = i (A^) 7475 (kjiiq) (H3) 



Both currents are invariant under global chiral transformations but they are not invariant under 
[/(l)-gauge transformations: 

^-n?^?, Tj^Tffl (H4) 

with some unitary function of discrete variables. In order to have gauge invariance we define a 
gauge field on the lattice 

m n = — ^4 — -IT H5 

1-It{(A )] + (A )]'} 

and 

Ujj> = — t-4 — -4 116 

" i-H(^i)?+w} 

where (Ai,iAq) are the two component electromagnetic vector potential, each of them satisfying 
the wave equation ( |22|) with M = namely, 

^VnAnVjAj - ^VjAjVnAr^j (Atf = (117) 

and similarly for {Aq)™. 

The gauge fields C/j 1 ™' are associated with the link between the points (j, n) — > (J, n') in the 
positive direction, and they transform under the gauge group as follows: 

u™' -» n] n u? n 'nf (us) 
up - <>]"r;>o;, (119) 

Inserting these fields in the vector and axial currents between fermion fields at separated points, 
we get 

l (—rn , n —rn n n +i , n +i -rre+1 ,-,-n+l.n ,n , T^i+l /n+l\ / 1on \ 



J4 



7 {^74^ n + ^74^. +1 ^ +1 + V7 t +1 74C/" +1 , i ^ + V? + i74^+i} (121) 



and similarly for and j\. Using ( |118| ) and (119) we can prove that all these expresions for the 
vector and axial currents are invariant under the gauge transformations ( |114| ). 

Now we want to calculate the vacuum expectation value of the divergence of the vector and 
axial current. We assume that this vacuum expectation value approches the corresponding non 
interacting fields for vanishing spacial separation of the fields ||. Thus we use solution of the free 
massless fermion field (^) with M = 0, namely, 

1 2 

1% = ~7= { U mCmf^j {k m ,UJ m ) + V m dbfj {k m , -0J m )j (122) 



'Ne 

m=-- 



with u r . 



and v r . 



mutation relations: 



Cm , C, 



and the operators c m and d m satisfying the usual anticom- 

= d mm ' (123) 



J mm' > 



rl tV 



Applying these operators to the vacuum of the Fock space, we get 

dtjo) = c m |0)=0for m>0 
d m \0) = cm\0) = Ofor m < 

Collecting these properties we obtain for the axial current 



1 A -5 1.* -5 



(124) 
(125) 



£ (0 1^7175^1 0) - 1 (0 |^ 7175 ^ 



+ 



+ 



4e 

4r 
1 

iVe 

1 

iVe 



n,n+l 



7475Vj+l 



0>C/ 



n+1 



N/2-1 

m=o \ 
JV/2-1 



1 - |to^ 



5i 



2 TU} m, 



1 + 2 e ^m 



m=0 



♦ t ( 

m=-7V/2 
\ -1 

+ £ 

m=-jV/2 



^7475^+1 



1 + §ro; m 
1 - |ro; m/ 



1 -f- 

1 ~2^km ^ 



0)C/^ +1 + c.c. 



'- / TT«,n+l Trn,n+l\ 

47 V i J 



4r 



-1 



-i-w/N 



+ 



2iV^f 1 - ^-A^^o) • - f A|(A ) J 2iV S enf 1 _ feA^A^ 
which in the limit, N — ► oo, e — > 0, r — > 0, becomes 

1 



t" _ i: 
'j 4 



v2n 
l) 2 j 



+ C.C 



dijl + <9 4 j 4 & 







- (diA) - 9oAi) = —F 10 

7T 7T 



(126) 



as required. Similarly we can prove 



1 A • 1 -A ■ 

-Ajji + -iA n j 4 
e r 











We come to the conclusion that our model leads to an interaction which is [/(l)-gauge invariant, 
but the divergence of the axial current gives in the continuous limit a photon mass, as expected by 
the axial anomaly. 

For the sake of simplicity in the notation we have worked all the formulas in (1 + l)-dimension. 
The expressions in this paper can be easily generalized to (3 + l)-dimensional lattice by the use of 
the promediate fields: 



V jl A jl Vj 2 A j2 Vj 3 A j3 (p(e 1 j 1 ,e 2 j 2 ,e 3 j 3 ,Tn) -> (f>(x 1 ,x 2 ,x 3 ,t) 



and similary for space an time differences. 



We have discussed symmetries and conservation laws for this scheme elsewhere] 10 1. 
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